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We present a new nucleon distribution amplitude which amalgamates features 
of the Chernyak-Ogloblin-Zhitnitsky model with those of the Gari-Stefanis model. 

(N 

\ This " heterotic" solution provides the possibility to have asymptotically a small ratio 

^ \ 

Q^l |GjJ^|/G^i < 0.1, while fulfilling most of the sum-rule requirements up to the third 

Q ' order. Using this nucleon distribution amplitude we calculate the electromagnetic 

and weak nucleon form factors, the transition form factor 7pA+ and the decay widths 

•T— I , 

^ ' of the charmonium levels ^Si, ^Pi, and into pp. The agreement with the available 

data is remarkable in all cases. 
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In recent years a number of authors |T|-^ have proposed various models for nucleon 
distribution amphtudes (NDA), based on hght-cone perturbative QCD in conjunction 
with QCD sum rules (SR) With varying degrees of conviction these models have been 

used in a series of analyses to calculate the electromagnetic and the weak p!0|JTl[ 

nucleon form factors, the transition form factor 7pA+ [P^^-[T3|,P|] , the cross section for proton 
Compton scattering at large momentum transfer and the exclusive pp decays of 



heavy quarkonia |1T|,|8|JT7[] . They may also be useful in determining the pion form factor at 
large spacelike via pion electroproduction |jl8 . 



Although these models incorporate essential ingredients for a unified description of per- 
turbative and nonperturbative aspects of the subnucleon structure, important questions still 
remain unanswered. One question focuses on the value of G^.j/G^j and the possibility that 
the electron-neutron differential cross section is dominated by G^, while G\j is asymptot- 
ically small or equivalently that ^ at all values [|19|. The Gari-Stefanis (GS) 



amplitude was constructed to account for this behavior of the form factors and gives good 
agreement with the latest high-Q^ SLAG data ||2^ at the expense that the moments (002) 
and (101) cannot match the SR requirements of Ghernyak and Zhitnitsky (GZ) in the 
allowed saturation range [^]. A second issue concerns whether higher terms in the Appell 
polynomial decomposition of the NDA ||5| are significant, a point that was raised in |^. A 
systematic investigation of this matter will be conducted elsewhere. In the present work we 
resort to second-order Appell polynomials. 

In Ghernyak, Ogloblin, and Zhitnitsky (GOZ) have recalculated the SR for the first- 
and second-order moments of the NDA and have derived for the first time SR for the third- 
order moments. Their new set of SR comprises 18 terms with restricted margins which 
comply with the results of the King-Sachrajda (KS) computation P] but disagree with those 



obtained on the lattice for the lowest two moments |22[. They have also proposed a new 
NDA which satisfies all, but 6 SR, whereas the GZ amplitude and the GS amplitude violate, 
respectively, 13 and 14 of the new SR. The KS amplitude provides almost the same quality 
as the GOZ amplitude with only 7 SR being broken. 



In the same authors pointed out that the GS model leads to a prediction for the 
^Si pp decay width which is about 50 times smaller than the experimental value. This, in 
connection with the strong violation of the SR for the moment (002) led Chernyak, Ogloblin, 
and Zhitnitsky to the conclusion that the GS model is unacceptable. On the other hand, 
they stressed [H,^ that "they have not succeeded in finding the model distribution amplitude 
which contains not higher than second-order Appell polynomials and gives |F"|/Ff < 0.4 
while fulfilling the SR." 

It is the purpose of this note to present a unique NDA which makes it possible to resolve 
all these problems. This novel solution turns out to be something of a hybrid between 
the COZ amplitude and the GS amplitude; thus the "heterotic" NDA. Specifically, we find 
I^mI/^m ^0.1 and a good agreement between the calculated proton form factor and the 
high-Q^ data pO|, whereas the fit to the SR has almost the same accuracy as the original 



COZ model. In addition, the heterotic distribution amplitude leads to predictions for the 
exclusive decays of the charmonium levels ^5*1, ^Pi, and into pp, which are in remarkable 
agreement with the existing experimental data. 

The theoretical basis for the description of hadronic exclusive processes within QCD 
is provided by the factorization property (see, e.g., [^), meaning that all soft (nonpertur- 
bative) effects can be absorbed into quark distribution amplitudes for the hadrons in the 
initial and final states, while the hard-scattering subprocesses can be calculated via pertur- 
bative QCD. The leading-order definition of the NDA, modulo logarithmic corrections due 
to renormalization, is 

$(x„Q2) = f\d'k^]ij{x„kf), (1) 

where the measure is [d^k^] = 167r35(2)(^3^^ k?)n-=i and ^pixi,kf) is the lowest- 

twist Fock-space projection amplitude for finding three valence quarks inside the nucleon, 
each carrying a fraction Xi = kf /p'^ (with p^ = p° ± p^) of the nucleon's longitudinal 
momentum p"*" and having relative transverse momentum 'k^^ . Although %p{xi,'k^^), and 
hence $(xj,Q^) cannot be calculated within perturbative QCD, a evolution equation 



can be derived from QCD perturbation theory ||^. Any solution of this equation can be 
expressed in the form 

°° (n (O'^W^"" ~ 

$(x„g2) = «l>..(x.)E^n(/x^)M^ $n(a:.), (2) 

n=0 \ Ctsl^/i ) ) 

where $as(xj) = 120x1X2X3 and {$„(xj)} are the eigenfunctions of the interaction kernel 
of the evolution equation, represented in terms of Appell polynomials (cf. Ref. The 
corresponding eigenvalues 7„ equal the anomalous dimensions of the lowest-twist three-quark 



operators with the appropriate baryonic quantum numbers ^ . 

The nonperturbative input enters Eq.(2) through the coefficients which represent 

matrix elements of appropriate three-quark operators (in the light-cone gauge = 0) 
interpolating between the proton and the vacuum: 

Here z is a lightlike vector with z'^ = 0, A^-^ is the proton spinor, and Jn denotes the "proton 
decay constant" . SR calculations make use of correlators between two of the 0(j;*i"2n3). 

(4) 

To determine the moments of the NDA , 

^(nin.na) ^ ^^^^ ^n, ^„3 ^ (5) 

Jo 

{[dx] = dxidx2dx36{l — Xi — X2 — X3)) a short- distance operator product expansion is per- 
formed at some spacelike momentum where quark-hadron duality is valid. By virtue of 
the orthogonality of the eigenfunctions the coefficients can be determined by 

inverting Eq. (5) upon imposing the SR constraints. The moments ^j^^^^^i^) terms of the 
coefficients Bn for n = 0, 1, ... 5 have been given in [^. Those for n = 6 . . .9 are calculated 
in E4 



Let us now outline our results. Treating the SR defined in (4) for ni + n2 + ris < 3 and 
m = 1 within the range estimated by COZ the coefficients i?„ for the heterotic solution 
are: Bq = 1, Bi = 3.4437, B2 = 1.5710, ^3 = 4.5937, B^ = 29.3125, and B5 = -0.1250; the 
value of Bq is due to the normalization of $Ar, i.e., /o^[(ia;] $Ar(xi) = 1. The exphcit form of 
K"' is 0: 

^n\x,) = ^as{xi){-19.773 + 32.756(xi - X3) + 26.569x2 

(6) 

+16.625x1X3 - 2.916x^ + 75.25x^} . 

The moments and the corresponding SR constraints are shown in Tab. |. There is a good 
overall consistency, with only 7 SR being broken. Given the fact that we have taken into 
account only the first 6 Appell polynomials to represent $Ar, the deviations are tolerable. 

The results for the magnetic form factor G'^j have been obtained using analytical ex- 
pressions given in and are plotted in Fig. [1| for the proton and Fig. ^ for the neu- 
tron. The data are from pD| , p5[ . The form-factor evolution with is due to the leading- 



order parametrization of the effective coupling constant ^^(Q^)- The evolution of the 
coefficients -B„ is a minor effect and has been neglected. Note that an average value 
tts(Q^) = [<ys{Q^ X 0.427)as(<5^ X 0.178)]^/^ has been used to account for the different 
virtualities of the involved propagators. Here and below two options are shown correspond- 
ing to two typical values of the scale parameter Aqcd, while the proton decay constant is 
taken to be |/^| = (5.0 ± 0.3) x IQ-^GeV^, as suggested by QCD-SR [|l],|,|. 

Using $5^f* we have calculated the electromagnetic A^ — A"^ transition form factor G\.^{Q'^) 
for two [0,0] recently proposed A"*" distribution amplitudes (Fig. 0) following ||12[. The solid 



and dashed curves refer to the model of Carlson and Poor (CP) ||2^ for Aqcd = lOOMeV" and 
Aqcd = 180Me\^, respectively. The dotted and dash-dotted curves are their counterparts 
for the model of Farrar et al. (FZOZ) In all cases the CP value |/a| = 11.5 x 

lO^^GeV^ has been used, which is within the spread of the FZOZ estimate. The data are 



compiled in ||28|. The evolution of G\^{Q'^) is governed by an effective coupling constant. 



which we take to be the average of the two couphng constants a^J^\Q'^) and a^^^(Q^) with 
arguments weighted by the virtuahties appropriate to each modeh a^'^^\Q'^) = x 
0.3773)as(Q2 x 0.1488)]i/2 and ap^^)(Q2) = [a,{Q^ x 0A6A3)as{Q^ x 0.1015)]i/2. if we 
take these predictions at face value, then the available data seems to favor the CP model 
(cf. [^). On the other hand, there is as yet no definite experimental evidence whether the 
Q^G\,j curve levels off or descends rapidly to zero, as predicted by COZ-type NDA |T^. A 
priori we have no reason to favor one option over the other (and in fact the recent analysis 
by Stoler |T5[ of the unpublished data of the SLAC experiment E133 points to the second 
possibility). Further exclusive experiments to measure G\.j at as high as possible are 
crucial. 

The calculation for the nucleon axial form factor QAiQ"^) according to yields at 
Q2 ^ lOGeV^, QV(Q^) = 0.90Ge1/^ for Aqcd = lOOMeV, and QV(Q^) = lAAGeV^ 
for Aqcd = ISOMeV. These results compare well with the value Q'^gA^Q'^) ~ 1.5GeV^ 
extrapolated from the data Also the ratio gA{Q^)/G\,j{Q^) ~ 1.19, in the region 

where the calculations can be trusted, is consistent with the (extrapolated) experimental 
value gA{,Q'^)/G\j{Q'^) ~ 1.35. As for the isoscalar nucleon form factor |Tl|, we find at 
g2 ^ IQGeV^, Q^G^a\Q^) = 0.83GeV^ for Aqcd = lOOMeV and Q^G^i'iQ^) = IMGeV^ 
for Aqcd = ISOMeV. Assuming isospin invariance, we combine these results with those for 
gA to obtain G^'^/G'^'* ^ 1.85, where G^'' is the isovector axial-vector nucleon form factor. 
If a dipole form G'x\Q^) = G^^\o) / {1 + M^s? , with G?(0) = 0.38 from SU{6), is 
used to describe the dependence of G^'' |]TI|], then, in the high-Q^ region, our model 
yields Mas = (1-15 - 1.22)GeV for Aqcd = lOOMeV and Mas = (1-27 - 1.37)GeV" for 
Aqcd = ISOMeV. These results might be relevant to studies concerning the strange-quark 
content of the nucleon. 

The last issue we address in this work are the exclusive decays of the charmonium 
levels ^Si, ^Pi, and into pp. Such calculations have been carried out by several au- 
thors [P,p|, p!7tpT] , p2[ within the QCD convolution framework. We here follow [§. We consider 
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first the two Xc(l-P) states. The branching ratio for the decay of the J^^ = l^"*" state into 
pp is given by 



BRl 



pp 



0.75 16n^ 


In 


ln(M/A) 729 


M2 



^3Pi ^ all J 

where M ^ 2mc 3GeV and A = OAGeV (the last value from |0-see also ||3 
nonperturbative content of Eq. (7) is due to and the amplitude for the process ^Pi 



(7) 
The 



Ml, which involves $Ar. Using (6) the calculation of Mi yields M^^* = 99849.6 and as a result 
BR{^Pi — > pp/^Pi — s> a//) = 0.77 x 10~^%, which is in accordance with the experimental 
value (0.5 - 1.0) x 10-^% |3|. 

The analogous expression to (7) for the J^*" = 2++ state has the form 

pp 



BRl 



0.85(7ra, 



16 


fN 


729 


M2 



2 ) 



(8) 



'2 (ill ^ 

which is Eq. (20) of with an obvious minor correction. For the heterotic NDA we find 
Ma = 515491.2. Setting a,(mc) = 0.210±0.028 (see third paper of |3|), we then obtain from 
(8) BR(^P2 pp/^P2 all) = 0.89 x 10~^% in remarkable agreement with the measured 
value (0.90 ^JJg ± 0.19) x 10^2% ^ . 

state into pp is 



The partial width of the J^*^ = 1" 

T{^Si pp) = {iras 



1280 \Uf 


fN 


2437r M 


M2 



M 
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(9) 



where determines the value of the ^S'l-state wave function at the origin. Its value 
can be extracted from the leptonic width T(^Si — > e+e^) = (4.72 ± 0.35)A;eV" ||35| via 
the Van Royen-Weisskopf formula. The result is |/^| = 383MeV with mj/^ equal to 
its experimental value. The heterotic solution leads to Mq = 13726.8. Then, using the 
previous parameters, it follows that r('^S'i — > pp) = 0.12keV. From experiment ||3^ 
it is known that T{pp)/Ttot = 2.16 ± 0.11 x lO'^ with Tm = (68 ± 10)keV, so that 
r(^S'i — > pp) = 0.15keV in excellent agreement with the model prediction. For the branching 
ratio we find BR{^Si pp/^Si all) = 1.76 x 10'^ or 1.40 x 10-^ if the new value 
Ttot = 85.5^ggA;eV^ is used. We emphasize that the model predictions for all considered 
charmonium decays are obtained with the same values of the various parameters. 



By considering a number of exclusive reactions involving the NDA $ jv , we have effected 
that a novel solution, which we call heterotic, leads to predictions which are corroborated 
by experiment. This solution is consistent with the SR requirements up to the third order 
and allows for the possibility to analyze the form-factor data with the assumption that 
asymptotically \G\j\/G\j < 0.1. We have pursued our approach, in spite the objections 
raised by Isgur and Llewellyn Smith and also by Radyushkin ||38|. We nevertheless 
believe that our model provides a useful and predictive tool for phenomenological studies. 
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FIGURES 

FIG. 1. The proton magnetic form factor calculated with the heterotic distribution amplitude 
in comparison with the data. 

FIG. 2. The neutron magnetic form factor calculated with the heterotic distribution amplitude 
in comparison with the data. 

FIG. 3. Comparison with available data of the transition form factor 7pA"'" calculated with 
the heterotic nucleon distribution amplitude and two different models for the A+ resonance, as 
explained in the text. 
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TABLES 

TABLE I. Moments ni + n2 + ria < 3 of the heterotic nucleon distribution amplitude $jv in 

comparison with the sum-rule constraints 



Moments {nin2ns) 


Sum 


rules 


^(nin2n3) 
^N/het 


(000) 




L 


1 


(100) 


0.54- 


-0.62 


0.572 


(010) 


0.18- 


-0.20 


0.184 


(001) 


0.20- 


-0.25 


0.244 


(200) 


0.32- 


-0.42 


0.338 


(020) 


0.065- 


-0.088 


0.066 


(002) 


0.09- 


-0.12 


0.170 


(110) 


0.08- 


-0.10 


0.139 


(101) 


0.09- 


-0.11 


0.096 


(Oil) 


-0.03- 


-0.03 


-0.021 


(300) 


0.21- 


-0.25 


0.21 


(030) 


0.028- 


-0.04 


0.039 


(003) 


0.048- 


-0.056 


0.139 


(210) 


0.041- 


-0.049 


0.079 


(201) 


0.044- 


-0.055 


0.049 


(120) 


0.027- 


-0.037 


0.050 


(102) 


0.037- 


-0.043 


0.037 


(021) 


-0.004- 


-0.007 


-0.023 


(012) 


-0.005- 


-0.008 


-0.007 
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